Sinusoidal excitations in reduced Maxwell-DufRng model 
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Sinusoidal wave solutions are obtained for reduced Maxwell-Duffing equations describing the wave 
propagation in a non-resonant atomic medium. These continuous wave excitations exist when the 
medium is initially polarized by an electric field. Other obtained solutions include both mono- 
frequency and cnoidal waves. 
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An atomic medium in general conditions is modelled by 
N-level atoms. In the two-level resonant approximation, 
the system is characterized by the Bloch equations, which 
is inaccurate in several physical situations [l], like 
dense atomic media and systems involving three or more 
level atoms. Thus the resonant model needs to be gen- 
eralized and extended to the non-resonant scenario. One 
of the well studied approaches is to consider the response 
of the medium as weakly nonlinear. Such situation leads 
to the Duffing oscillator model, where the nonlinear re- 
sponse of the medium is assumed to be cubic. This is the 
simplest generalization of the Lorentz model. 

On the other hand. Maxwell wave equation, for a lin- 
early polarized light, allows propagation in both the di- 
rections. However, this can be approximated to unidirec- 
tional wave propagation when anharmonic contribution 
to the polarization is very small. As a result, the wave 
equation reduces from second order to a first order equa- 
tion. For a non-resonant medium, this approximation 
results in the reduced Maxwell-Duffing model (RMD). 

Different excitations in non-resonant atomic media are 
currently attracting considerable attention, because of 
their relevance to ultra-short regime. Detailed reviews of 
various aspects of non-linear excitations in atomic media 
can be found in [1, In case of two level atoms, the lo- 
calized soliton solutions of the Maxwell-Bloch equations 
explained the physical phenomenon of self-induced trans- 
parency 0. In the same system, general cnoidal waves 
have also been found as exact solution [1, 0, [1] ■ It was 
observed [§] that, these waves can be naturally excited in 
the presence of relaxation. Such shape preserving Jacobi 
elliptic pulse train solutions have been experimentally ob- 
served [lol . More general periodic solutions in multi- level 
systems have also been reported 

In case of Maxwell-Duffing model, a class of exact lo- 
calized soliton solutions have been recently obtained ■ 
We present here mono frequency, sinusoidal wave exci- 
tations for RMD system. This excitation exists only in 
the presence of a polarizing background. General cnoidal 
wave solutions are found both with and without back- 
ground. 



Below a brief summery of the reduced Maxwell-Duffing 
model is presented after which a procedure to find ex- 
act solutions of this system through a fractional linear 
transform is outlined. We then present the novel si- 
nusoidal wave solutions including single frequency and 
cnoidal waves. 

Reduced Maxwell- Duffing Model: 

The propagation of electromagnetic waves in a medium 
is described by the wave equation: 



(1) 



where P is polarization of the medium. For unidirec- 
tional wave propagation, the above equation can be re- 
duced to a first order equation. 
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In the Duffing oscillators model, the nonlinear response 
of the medium is cubic. The corresponding equation for 
the motion of electrons in the presence of an external 
electric field is given by 
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Here X represents the displacement of an electron from 
its equilibrium position, ojq is the oscillator frequency, K3 
is anharmonicity constant, and m is the effective mass 
of the electron of charge e. The medium polarization is 
defined as P = neX, where n is the number density of 
the oscillators in the medium. 

We choose new variables r = z/l, x = LOQ(t — z/c) and 
normalize the independent variables as 

~e = E/Ao,q = X/Xo. (4) 
In terms of the new variables, Eq. ([3]) then takes the form 



e, 



(5) 



where 2fi 



K3X^/0l and Aq 



mujQXo/e 



'^{2fi/\K3\y/^ . Xq can be expressed as Xq 
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(2/xa;o/|K3|)^/^. Similarly Eq. ^ is transformed to 

de dq 
dr dx 



(6) 
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FIG. 1: Mono-frequency solutions for /i = 0.5 with (a) a = 1.6 and (b) a — 3.0. The second case imphes A — 0. 



Here I is defined as 

= 27rneV("ic0o) = 0p/2c0o, (7) 

where LUp = (47rne^/m)^/^ is the plasma frequency. 
Eq. (O and ([6]) together are caUed reduced Maxwell- 
Duffing equations. For finding propagating solutions, one 
defines 



Fractional Linear Transform and the solutions: 

The above Eq. pop is in the form, which can be ob- 
tained from the real part of the non-linear Schrodinger 
equation with a source. For finding out explicit solutions 
we consider Eq. (fTU|) : 

q" + gq^ + eq^C, (11) 



ri — X — T / a, 



(8) 



where a is related to the velocity of the pulse. Eq. ([6]) 
can be integrated with respect to the single variable rj to 
yield 



e(t, x) = aq{t, x) + C, 



(9) 



where C is a constant, which signifies the background 
electric field, when electron amplitude q goes to zero. 
Non-linear equation of motion in Eq. ([5]) then takes the 
form: 



^ + {l-a)q + 2fiq^ = C. 



(10) 



For a wave propagating in the right direction, a > 1 and 
fi > 0, which is considered below. The other propagation 
direction can be like wise studied. 



provided g = 2fi and e = (1 — a). Prime indicates differ- 
entiation with respect to rj. It is known earlier ^13] , that 
this equation can be connected to the elliptic equation 
/" + af + bf^ = through the following fractional linear 
transformation (FT): 



qiv) 



A + Bf{r],mY 



(12) 



l + Df{Tj,m)\ 
where A, B and D are real constants, S is an integer, and 
/(?7, m) is a Jacobi elliptic function, with the modulus 
parameter m (0 < m < 1). It can be shown, that 6 = 2 
is the maximum allowed value. Here we concentrate on 
the periodic solution of Eq. pT|) for (5 = 1 and q{r], m) = 
cn{ri,m). Solutions for other elliptic functions also can 
be studied in a similar way. 

/. General solution: 

i) The consistency conditions for m = are given by 



A^g + 2D{AD - B) + Ae - C = 0, (13) 

SA^Bg + AD{1 + 2e) + B{e - 1) - 3CD = 0, (14) 

SAB^g + AD^{e - 1) + BD{1 + 2e) - 3CD^ = 0, (15) 

B^g + BD^e - CD^ = 0. (16) 

I 



It should be pointed out that cn(?7, 0) = cos{ri). One implies only a constant or trivial solution and is not con- 
can see from the above FT (Eq. ([12])) that AD = B sidered here. An {AD — B) factor can be taken out of 
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FIG. 2; Cnoidal wave solution in presence of source, (a) (Dotted line) g — 2.0, (solid line) g — 7.0 and (dashed line) g = 30 for 
e = -5.0; b) (Dotted line) e = -0.001, (solid line) e = -2.0 and (dashed line) e = -5.0 for 3 = 5. 



all the conditions by tactically using the first consistency 
in Eq. (|13p . Other conditions were used to solve for the 
unknowns A, B and D. The source term (C) can be 
determined from the first condition itself. Thus, the so- 
lution is expressed as 



A + Bcos{r]) 



1 + Dcos{rf) 
where, A^ B and D respectively are 



(17) 



and 




-(l + 2e) 



6g 



(18) 



After solving the solution parameters, the source term or 
the constant electric field part can be determined from 
Eq. HT 



C 



(1 + 2e) / (1-e) 
3 V 35 



(19) 



As has been mentioned earlier, for wave propagation in 
the right direction, /i > and a > 1, implying g > and 
e < respectively. The parameter values exhibited in 
Eq. yield the domain of the solutions: g > Q and e < 
— 1/2, where A, B, D are real and C is a positive quantity. 
It is worth pointing out that all the solutions are non- 
singular in nature. This is because the magnitude of D 
is less than unity. The solutions are depicted in Fig.[T](a) 
for fi — 0.5 with a = 1.6. The dotted Hne corresponds to 
the solution for positive value of D and the solid line is 
for the negative one. The first plot is with a background, 
i.e., A^O. 

ii) The consistency conditions, (jlSlllGp support solu- 
tion for A = if the source is non-zero. In this case, 
e = — 2 and the solution is written as 



^ 1 / cos{t]) 



y/g V\/2± cos(?7) 



(20) 



which is plotted in Fig. [T] (b) for /i = 0.5 with a — 3.0. 

iii) Although a wide class of localized pulse propaga- 
tion is analyzed in for different initial boundary con- 
ditions, for completeness we indicate the same for m = 1 
in the fractional transform. This is expressed as 



A + Bsech{r]) 
1 + Dsech{ri) ' 



(21) 



where 



A = ± 



3g 



'(l-2e)' 



B^± 



2(e-2)2 
3s(l-2e) 



with = -(1 + e)(l - 2e)2/(27.g). This solution exists 
for 5 > and e < — 1. For e — 2, the solution goes to 
the one, with B = Q. All the solutions are non-singular 
except for B — Q and D — — v^, which implies a singular 
one, signifying self focussing effect. 

iv) We now analyze the nature of the solutions in the 
absence of a polarizing background: C = 0. The resul- 
tant dynamical equation is the real part of the well known 
non-linear Schrodinger equation. In this case, the peri- 
odic solution for the electron a mplitu de is of quadratic 
fractional type [5 = 2): A = ^J-^/g, B ^ 0, D ^ -2 
and e — 4. The solution is singular one, implying an in- 
stability in the electron amplitude or the self focussing 
of the electric field. This happens for any value of /i < 
and for a = —3. 

v) In addition to the above mentioned mono-frequency 
solutions, there are other conoidal wave solutions. As an 
example, when m = 1/2 the consistency conditions are 

A^g + D{AD - B)+Ae-C ^0, 
3A^Bg + {2AD + B)e - 3CD = 0, 
SAB'^g + D{AD + 2B)e - iCD"^ ^ 0, 
B^g + BD^e + {AD - B) - CD^ = 0. 
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These equations lead to, 




where E ^% + 2e^ - 3\/9 + 4e2 and F = V^Ti? - 3. 
E and F are positive for all e ^ 0. These conoidal waves 
exist in the domain e < and g > 0. Solutions are 
picturised in Fig. [21 where the variations of the electron 
amplitude with g and e are displayed. As /i increases the 
amplitudes diminish for a fixed value of a = 6 (Fig.[2fa)). 
Fig. ^h) shows the nature of the solutions with a for 
certain value of fi = 2.5. 



In conclusion, novel periodic solutions have been ob- 
tained for the reduced Maxwell-Duffing model, which are 
unidirectional, sinusoidal and of mono-frequency charac- 
ter. These solutions complement the localized soliton 
solutions known in the literature and occur only when 
a back ground electric field is present. General cnoidal 
wave solutions are found both with and without back 
ground fields. It will be interesting to study the nature 
of the waves when Maxwell equation is modified by a 
non-local dispersion term [l^ . Polarization properties 
of the propagating modes should also be investigated in 
this non- linear medium p^ . as also the nature of the 
waves when slowly varying amplitude approximation is 
not valid p^ . 
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